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Annihilation of clumped superheavy dark matter provides an interesting explanation for the origin of ultrahigh 
energy cosmic rays. 

The predicted anisotropy signal provides a unique signature for this scenario. 



1. INTRODUCTION 

Explaining the origin of cosmic rays at ener- 
gies above the Greisen-Zatsepin-Kuzmin bound 
Eqzk — 4 x 10 19 eV is an interesting and puzzling 
scientific challenge for several reasons: 

AGASA observes a flux of order of one cosmic 
ray event per 100 km 2 per year at energies above 
10 20 eV [Q , thus confirming the flux seen by previ- 
ous smaller ground based arrays on the northern 
hemisphere]], see e.g. ||. 

On the other hand, energy loss of electromag- 
netically interacting particles due to scattering 
with the CMB background indicates that these 
ultrahigh energy cosmic rays (UHECRs) should 
originate within a radius of order 100 Mpc j| [|, 
and sufficiently powerful astrophysical accelera- 
tion mechanisms are not easy to find. 

This has motivated proposals that superheavy 
dark matter may be the origin of ultrahigh energy 



1 But see |g| for a recent discussion of the pertinent prob- 
lem to reconstruct primary cosmic ray energies from air 
showers. A new analysis presented there lowered the 
estimate on the energy of the highest energy event ob- 
served with the Haverah Park Array from above 10 2t) eV 
to 8.3 X 10 19 eV and also indicates a smaller flux at the 
high energy end of the spectrum. Of course, the very ex- 
istence of UHECRs is not questioned by these findings. 



cosmic rays: 

This includes in particular Hill's proposal of 
monopolonium decay , the decay of superheavy 
dark matter relics proposed by Berezinsky et al. 
^ and by Kuzmin and Rubakov [^) , and most re- 
cently our proposal of collisional annihilation of 
superheavy dark matter ]lO| , pd| . While decay sce- 
narios for UHECRs should imply domination of 
the galactic center in the distribution of arrival 
directions, the proposal of collisional annihilation 
predicts that with sufficient statistics from the 
upcoming Auger observatory a clumpy structure 
should emerge in the galactic source distribution: 
Unitarity bounds on reaction cross sections im- 
ply that a smooth superheavy dark matter dis- 
tribution in the galactic halo would not yield the 
observed flux through collisional annihilation, but 
the overdensity in dark matter clumps in the halo 
may well account for the observed flux Jl(| . 

The common virtues of these top-down scenar- 
ios are: 

• a natural explanation of the seemingly random 
distribution of arrival directions; 

• avoidance of collisional damping and the GZK 
bound due to origination in the galactic halo; 

• no need for an extremely powerful astrophysical 
acceleration mechanism. 



2 



Very high-energetic neutrinos both from super- 
heavy dark matter decay and from annihilation 
have very recently been discussed in p2[ , and 
some recent overviews of top-down scenarios for 
UHE cosmic rays can be found in [jl3| — |l5|| . 

Superheavy dark matter particles cannot arise 
as a consequence of standard thermal freeze out, 
but it has been pointed out that dark particles 
might have never experienced local chemical equi- 
librium during the evolution of the universe, and 
that their mass may be in the range 10 12 to 
10 19 GeV, much larger than the mass of ther- 
mal wimps p6[-pj|. Since these wimps would be 
much more massive than thermal wimps, such 
superheavy dark particles have been called WIM- 

PZILLAS 

wimpzillas may be created during bubble col- 
lisions if inflation is completed through a first- 
order phase transition [^0 21|; at the preheat- 
ing stage after the end of inflation with masses 
easily up to the GUT scale [^2| or even up to 
the Planck scale ^,^4|; or during the reheat- 
ing stage after inflation [Q. WIMPZILLAS may 
also be generated in the transition between an in- 
flationary and a matter-dominated (or radiation- 
dominated) universe due to the "nonadiabatic" 
expansion of the background spacetime acting on 
the vacuum quantum fluctuations [116 25 26J . The 



distinguishing feature of this mechanism is the ca- 
pability of generating particles with mass of the 
order of the inflaton mass even when the parti- 
cles only interact extremely weakly (or not at all) 
with other particles. 

The mass m of unstable dark matter which 
decays through gravitational couplings is con- 
strained by the fact that the lifetime is only of 
order r < m^ lanck /m 3 |2^^], and therefore un- 
stable heavy relics must decay through strongly 
suppressed channels which effectively correspond 
to sub-gravitational couplings. This provided one 
motivation for the collisional annihilation pro- 
posal in [|l0]], since the need for 2-particle colli- 
sions implies a larger lifetime of annihilating wim- 
pzillas as opposed to decaying superheavy par- 
ticles. 



2. THE FLUX FROM ANNIHILATION 
IN DARK MATTER CLUMPS 

The spectral fluxes from decay or annihilation 
of WIMPZILLAS of mass Mx and particle density 
nx(r) are 



and 
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respectively. The shape of the spectrum is de- 
termined by the function dM(E, Ej et )/dE, which 
gives the number of sufficiently long-lived parti- 
cles per energy interval emerging from a primary 
jet of energy E^ . Birkel and Sarkar and Berezin- 
sky and Kachelriess have used Monte Carlo sim- 
ulations to calculate the spectrum from decay of 
superheavy relics |2^|3(| . Fodor and Katz as well 
as Sarkar and Toldra employed numerical integra- 
tions of the DGLAP evolution equations pl[p2| . 
The resulting spectrum in the interesting range 
between 10 19 eV and 10 20 eV is similar to the 
spectrum from the old modified leading log ap- 
proximation (MLLA) |53j, which was employed 
by Berezinsky et al. in their initial proposal of 
UHE cosmic rays from decay of superheavy dark 
matter decay ||. 

We had also used the MLLA limiting spectrum 
in our proposal of wimpzilla annihilation 
(Fig-i- 

The annihilation cross sections from the nor- 
malization in Fig. [l] were calculated for a smooth 
NFW |Q background halo. These cross sec- 
tions would be in conflict with unitarity bounds 
on reaction cross sections for particles of mass 
Mx > 10 12 GeV and therefore the flux from wim- 
pzilla annihilation in the smooth halo compo- 
nent will (most likely) be negligible^. However, a 
fraction around a few percent of the dark matter 



2 There is also the possibility that WIMPZILLAS form meta- 
stable bound states, similar to monopolonium. Then the 
contribution from the smooth halo component could be 
important and would resemble the anisotropy signature 
expected from monopolonium or superheavy dark matter 
decay. 
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Figure 1. UHE cosmic ray spectra from su- 
perheavy particle annihilation (upper panel) or 
decay (lower panel). 



in the halo is expected to exist in dense isother- 
mal dark matter sub-clumps^], and due to large 
overdensities in the cores wimpzilla annihila- 
tion in these sub-clumps can comfortably gen- 
erate the observed UHE flux with annihilation 
cross sections well below the unitarity bounds 
|i"o[ . This does not change the good fit of the 
annihilation model to the spectrum, which is a 
concequence of dN(E , Mx) / dE , since the differ- 
ence between events originating in the smooth 
and clumpy components affects only the overall 
factor Jd 3 rnx{r) 2 \r Q -r|~ 2 in (||). 

Due to the negligibility of a possible contribu- 
tion from the smooth halo component a clumpy 
source structure will be an unmistakable signa- 
ture of collisional wimpzilla annihilation as a 

3 See e.g. js 5|] for evidence for sub-clumps from TV-body 
simulations, and pj) for evidence from gravitational lens- 
ing. 



source of UHE cosmic rays. In particular, nearby 
sub-clumps will appear as hot spots in the sky 
(see Fig. 4 in Ref. go|). 

To calculate the expected overall anisotropy 
(after averaging over the clumps in the case of 
annihilation) one can integrate the local decay or 
annihilation rates entering in (|l|) and (||) along 
lines through the galactic halo, taking into ac- 
count the galactic halo density and the expected 
sub-halo distribution, respectively. 

The left panel in Fig. 2 shows the result for 
annihilation in a relatively broad subhalo distri- 
bution with the sub-clump density dropping off at 
a characteristic radius Rf = 20 kpc (see for a 
discussion of the expected sub-halo distribution) . 





Figure 2. The left panel shows the envelope of the 
UHE flux from annihilation in a galactic subhalo 
distribution as a function of the angle towards the 
galactic center, whereas the right panel shows the 
UHE flux from decay in a galactic NFW profile. 



3. CONCLUSION 

All top-down scenarios for generation of UHE- 
CRs from superheavy dark matter predict a 
large photonic component in cosmic ray primaries 
above Eqzk — 4 x IO 19 eV, and a more or less 
dominant role of the galactic center: 

Both annihilation in bound states (monopolo- 
nium or " wimpzillonium" decay) and decay 
of superheavy relics imply that the dominant 
smooth dark matter halo distribution would be 
the primary source of UHECRs, whence the dis- 
tribution of UHECR arrival directions should 
track the overall dark matter halo. 

Collisional annihilation, on the other hand, 
would almost exclusively appear in the cores of 
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local dark matter overdensities, and therefore 



would primarily track local overdensities. Here 
some domination of arrival directions pointing 
back towards the galactic center is expected be- 
cause the sub-clump population is expected to 
drop off with distance from the galactic center. 
A nearby clump out of the line of sight to the 
galactic center could dominate the local signal 
and appear as a hot spot. 

The southern site of the Auger observatory sees 
the galactic center, and therefore the prediction 
of a clumpy source distribution with some accu- 
mulation towards the galactic center should be 
testable after about three years of full operation 
of the southern site, which should yield of the or- 
der of 100 UHECR events with energies at and 
beyond 10 20 eV. 
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